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Abstract. In this paper we present a brief description of a ladder operator formalism 
applied to orthogonal polynomials with discontinuous weights. The two coefficient 
functions, A n (z) and B n (z), appearing in the ladder operators satisfy the two 
fundamental compatibility conditions previously derived for smooth weights. If the 
weight is a product of an absolutely continuous reference weight wq and a standard 
jump function, then A n {z) and B n (z) have apparent simple poles at these jumps. We 
exemplify the approach by taking wq to be the Hermite weight. For this simpler case 
we derive, without using the compatibility conditions, a pair of difference equations 
satisfied by the diagonal and off-diagonal recurrence coefficients for a fixed location 
of the jump. We also derive a pair of Toda evolution equations for the recurrence 
coefficients which, when combined with the difference equations, yields a particular 
Painleve IV. 
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1. Introduction 

We begin this section by fixing the notation. Let P n (x) be monic polynomials of degree 
n in x and orthogonal, with respect to a weight, w(x), x G [a, b], 
b 

P m {x)P n (x)w(x)dx = h n 5 mtn , (1) 

where h n is the square of the weighted L 2 norm of P n . Also, 

P n (z) = z n + p 1 (n)z n - 1 + ... . (2) 

For convenience we set w(a) = w(b) = 0. The recurrence relation follows from the 
orthogonality condition: 

zP n (z) = P n+1 (z) + a n P n (z) + /3 n P n -i(z) n = 0, 1, ... , (3) 

where (3 P^i(z) := 0, the a n , n = 0, 1, 2, . . . are real and (3 n , n = 1, 2, . . . are strictly 
positive. 
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In this paper we describe a formalism which will facilitate the determination of the 
recurrence coefficients for polynomials with singular weights. Two points of view lead 
to this problem: On one hand the X-ray problem pQ of condensed matter theory, on the 
other hand related problems in random matrix theory which involve the asymptotics of 
the Fredholm determinant of finite convolution operators with discontinuous symbols [Bj. 
This paper is the first in a series that systematically study orthogonal polynomial where 
the otherwise smooth weights have been singularly deformed. The ultimate aim is the 
computation for large n of the determinant D n of the n x n moments or Hankel matrix 
(t JL j+k)'jfc=o with moments fij := f^w(t)Pdt where j = 0, 1, . . ., thereby doing what has 
been done for the determinants of n x n Toeplitz matrices with singular generating 
functions [2J. 

The deformed weight with one jump is 

w{x) = w Q {x)(l-^ + P6{x-x)\ , (4) 

where a < x < b is the position of the jump, 9(.) is the Heaviside step function and the 
real j3 parametrises the height of the jump. More generally, we take wj > to be the 
canonical jump function wj(x) = 1 + J2j=i Aj9(x — Xj) and w(x) := Wq(x)wj(x). 
The actions of the ladder operators on P n (z) and P n _\(z) are 

+ 5n(Z) ) Pni - Z) = PnMz)Pn~l{z) (5) 

(J- - B n (z) - y/ (z)J P n ^{z) = - A n ^(z)P n (z) (6) 

Mz) :=Y. — + T f y ' o{z) - y ' o{y) ^(yMy)dy (7) 

j=i n j a ^ y 

B n (z) := £ -^4- + j±- f ^ z) -^ y) P n .MPn{y)w{y)dy (8) 

j^Z-Xj hn-i J a Z-y 



where Wq(x) =: exp(— v (x)). If w(a) and w(b) are non-vanishing one must add 



P n (y)P n ^(y)w(y) 



y=a 



and 

K-i\v - z) 



y=b 



y=a 



h n (y - z) 

to (JJJ) and (jHJ) respectively. 

Now A n (z) and B n (z), the coefficient functions in the ladder operators, satisfy 
identities analogous to those found for smooth weights E] : 

B n+1 (z) + B n {z) = (z-a n )An(z) -v' Q (z) (9) 

B n+1 (z) - B n (z) = ^j±Wj zMnzM ~ 1 (io) 

Z Oi n 

R n (xj) = -R n (Ai, Ax, Xj) := jf-P^x^woixj) (11) 
r n (xj) = r n (Ai, A Af ,x j ) := — ^-P n (%)P n _i(%)w (%) (12) 



Orthogonal polynomials with discontinuous weights 



3 



The derivation of will be published in a forthcoming paper where the weight 

has several jumps and w (x) = (1 — x) a (l + x) 13 , x 6 [— 1, 1] is the Jacobi weight. 

Multiplying the recurrence relation (jHJ) evaluated at z = Xj by AjWo(xj)P n (xj) and 
noting (jllj) as well as (fl"2"j) we arrive at the universal equality 

(xj - a n )R n (xj) = r n+1 (xj) + r n (xj) . (13) 

Similarly, squaring r n (xj) we find a second universal equation 

r 2 n {xj) = /3 n i2 n (%)i? n _i(%) . (14) 

Note that in the expressions for A n (z) and B n (z) only v , the "potential" associated 
with the smooth reference weight, appears. The discontinuities give rise to R n (xj) and 

r n{Xj). 

It is clear from (JJJ) and (JBJ) that if v' (z) is rational, then A n (z) and B n (z) are also 
rational. This is particularly useful for our purpose which is the determination of the 
recurrence coefficients, for in this situation by comparing residues on both sides of (JHj) 
and (fTUj) we should find the required difference equations 

In the following section the above approach is exemplified by the Her mite weight, 
wq{x) = exp(— x 2 ) and w given by (J3J. It turns out that in this situation a n and (3 n are 
related to R n and r n in a very simple way. 



2. Hermite weight with one jump 

Now, wo(x) = exp(— x 2 ), so that v (x) = x 2 , and w(x) as in (jlj). Also, 

R n := R n {x) = (3S^w (x) (15) 

r n := r n (x) = (J — w {x) , (16) 

which are independent of the particular choice of Wq and 

A n (z) = + 2 and B n (z) = , (17) 

z — x z — X 

particular to v (x) = x 2 . Note that P, } - X ^ IS the value of Pj (z, x) at z = x. 

Instead of proceeding with the full machinery of and ()10|) we take advantage of 

the fact that Mq{x) = x 2 . From orthogonality and the recurrence relation, we have 

I f°° ^ R n (x) 

a n (x)h n (x) = - v' (x)P 2 (x,x)w(x)dx = — ^ — h n (x) (18) 



2 „ 

by integration by parts. The string equation, 

POO 

nh n -i(x) = / P' n {x, x)P n -i(x, x)w(x)dx , (19) 



is an immediate consequence of the orthogonality condition. Again, an integration by 
parts and noting that (3 n = h n /h n _i produces 

PnW = \ + ^ " ( 2 °) 
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It should be pointed out here that in general neither the string equation (|19|) nor (|18j) 
will provide the complete set of difference equations for the recurrence coefficients which 
can be seen if wq were the Jacobi weight. In such a situation the compatibility conditions 
(jnj) and (JTDJ) must be used. 

Now (H3J) and (JH|) become 

- = a n {x - a n ) . (21) 

and 

r 2 n = 2(n + r n )a n a n -i . (22) 
Equations (|21|) and ()22|) . supplemented by the initial conditions 



a (x) = ^exp(-x 2 ) 



POO 

[1-/3/2)^ + (3 / exp(-t 2 )dt 

J X 



and r (x) = 



can be iterated to determine the recurrence coefficients numerically. Also, explicit 
solutions to (}2Tj) and (}2"2"j) can be produced for small n. 

3. Derivative with respect to x and Painleve IV 



If (121 1) and ()22)1 are combined with the evolution equations to be derived in this section, 
the Painleve IV mentioned in the abstract is found. We begin with the L 2 norm h n (x), 
Equation (JTJ), which entails 

POO 

d x h n = -(3 I P 2 (x)w (x)5(x - x)dx = -f3P 2 (x)w (x) = -h n R n = -2h n a n . (23) 



and thus d$ In (3 n = — 2(a n — a n _i) since In (3 n = hih n — In h n _i. With (j2T!|) . 

-d s r n = a n _i - a n (24) 



2(n + r n 

which is the first Toda equation. Taking the derivative with respect to x of (|T| 
at m = n — 1 and using the definition of the monic polynomials (j2J) then gives 
dxPi(n, x) = r n since Pi(n, x) — pi(n + l,x) = a n (x) is an immediate consequence 
of the recurrence relation. Therefore 

d x a n = r n - r n+1 , (25) 

the second Toda equation. 

Eliminating r n+1 from ()21|) and the second Toda equation (|23j) . gives r n in terms 
of a n and 9^a n : 

r n = a n {x - a„) + \d~ x a n . (26) 

Using the first Toda equation (j2^|) to express a n _i in terms of a n and (9 £ r n and 
substituting (|2^j) into (|2^j) produces a particular Painleve IV jSj, 

= + - 8xa 2 n + 2(x 2 - (2n + l))a n , (27) 

2a n 
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which can be brought into the canonical form with the replacements a n — > a n /2 and 
x — > —x. Equation (|27jl is supplemented by the boundary conditions limg^-too a n (x) = 
0. In a recent paper [HJ, a Painleve IV was derived for the discontinuous Hermite weight 
using an entirely different method. 

Based on (J2H|) and (f23|) the derivative of the logarithm of the Hankel determinant 
D n = YYj=o hj can be computed as 

n— 1 n— 1 n— 1 p"2(~ ~\ 

d x \nD n = -2^a J = -^i? J = -/jexp(-x 2 )^ (28) 

<9 2 lnP n = 2r n . (29) 

where (|15)) has been used in the first line, which can be summed by the Christoffel- 
Darboux formula, 



P j (x,x)P j (y,x) _ P n (x,x)P n - 1 (y,x) - P n (y,x)P n _i(x,x) 



n-l 

y 

p£ hj{x) (x-y)h n „i(x) 



In the limit x — > y, we find, in general, 



/i n _i(x) 



1 P?(x, x) 



= /3 n A n (x)P^_ 1 (x,x) + A n _ 1 (x)P%(x,x) - (2B n (x)+y' (x))P n (x,x)P n _ 1 (x 1 x) 
using the the ladder operators (0) and (JBJ). With (JTTj) this entails 

n ~ l pi( x x) 1 

V \ \:, J = t pr (2/5„P n 2 „!(x, x) + 2P 2 (x, x) - 2xP n (x, x)P n _i(x, x)) 

hj{x) K-i{x) 

+7 FT (^fin-P„ 2 -i(i, ^) + P„-iP 2 (x, x) - 2r n P n (x, x)P„_i(x, x)) 

(X - X)hn-i{x) 

The apparent pole at x = x can be shown to have vanishing residue by considering 
/3exp(— x 2 ) x (residue): 

/?exp(-x 2 )^„ J R n P 2 „ 1 (x,x) + P n _iP 2 (x,x) - 2r n P n (x,x)P n _i(x,x)j//i n _i 
= 2(n + r n )a n a n -i + 2(n + r„)a n _ia n - 2r 2 = 2(2(n + r„)a„a„_i - r 2 ) = , 

where the last equality is due to (|22jl. A further regular term can be found as a 
contribution from the Taylor series of P, (x, x) about x = x namely, 

2-^P n P n _ 1 (x, x)P,U(x) + 2^P n (x, x)P' n {x) 

-2-^d x (P n ^(x,x)P n (x,x))\ x= - (30) 
ft-n-1 

where Pj(x) := d x Pj(x, x)\ x=x . Using the fact that f3 n = h n /h n -\ and (|T5j) . we see that 
the first two terms of (}30j) combined into 2(r n /h n -i)d x (P n (x,x)P n -i(x, x))\ x=x cancel 
the third. We are therefore left with the regular term: 

d x lnP n = -fi( 2/? /"T l(M) + 2 SMS. _ 2 - iVM)Pn-i(s,sn exp( _ £2) 

V "7i-l "n-l "n-l / 

= 2xr n - 2(n + r n )(a n + a„_i) . (31) 
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Using (El-flU it follows from 2xr' n - 2r' n (a n + a n _i) - 2(n + r n )(a' n + a' n _ x ) = 
A(n + r n ) (r„_i — r n — o^ n _i) = that (}3~T]) reproduces the second derivative (|2^|) . 
dl\nD n = 2r n . 

Let F n (x) := — \a.D n [x) be the free energy. Expressing (}3*T]) in terms of et n (x) and 
finally in terms of fy n (x), where a n (x) =: |(\l/„(x)) 2 , the free energy reads 

F n (x) - F n (-oo) = f I ((An + 1 - 2x 2 )^ 2 n (x) + ?>(3x^ n (x) - (3 2 ^ n (x)) dx 

2 J -OO 

+ f*n(5)<(*) , 

where F n (— oo) is the free energy corresponding to w(x) = exp(— x 2 ). Note that 
exp(— F n (±oo)) = 2ix n l 2 n)b=i ^2^! which gives rise to the sum rule 

! ((An + 1 - 2x 2 )ty 2 n (x) + 3(3x^ 4 n (x) - (3 2 ^ n (x)\dx = . 

J —oo 

With a minor change of variables (j2Uj) becomes the Toda molecule equation. First 
we note that 

r n (x) = 2/3„(x) - n and /3 n (x) = — = — . 

K^x) (D n (x)Y 

Defining D n (x) =: exp(— nx 2 )D n (x) it then follows 

a, 2 lnD ra = 4 ^ +A - 1 . 

We may express a„(i) in terms of the derivatives of the free energy, by noting that 
—F' n = 2xr n - 2(n + r n )(a n + a n - X ) and -F% = 2r n . One finds with (p4]l 

/~\ _ jjt ~ ^-^n ~l~ ^jn _ ^ _ \ fn ^/n 

4n _ 2 ^ "2 4 /« 
where F n (5) =: / n (x) + ni 2 . 

4. Asymptotics and Numerics 

For x = we find the asymptotic expansion 

a n (0) = -ifl + (-l) n sin(261nn + B) (32) 



2n 

- — (1 + (-l) n (sin(2Mnn + B) - 46cos(261nn + B))) + 0(n~ 2 ) 

b 2 / \ 

r n (0) = -&(-l) n cos(26mn + 5) 1 + sin 2 (261nn + B) + 0(n- 2 ) (33) 

2n V / 

guided by the numerics on the difference equations (jSj) and (j22J). The constant 6 is 

given by 

»:= i-lnfil^ with fe(-l,l) 



2vr \1-P/2J 2 




Figure 1. Behaviour of 6 _1 (n/2) 1 / 2 a n for (3 = 3/2 and various (small) values of x 
from a numerical iteration of (|21l) and (|22|l . The top panel shows additionally Equation 
(13211 as a solid line passing through the open circles. The other two panels show the 
approximation of the rescaled a n {x) using the first three (two) terms of a Taylor series 
around x — as dashed (dotted) line. The behaviour of r„ looks likewise. The 
numerical data is pruned differently in the three cases to avoid artefacts. 

and B is a phase independent of n. Unfortunately, the formalism developed in this 
paper does not seem to shed any light on its delicate dependence on (3. 

It can be verified by a direct calculation that (jH2j) and (jHHJ) satisfy (J2H) and (j22J) 
to order 1/n 2 . The top panel of figure ^ shows a comparison between the numerical 
results and the above asymptotes for suitably rescaled a n . In principle, a n and r n can 
be determined analytically to any order in n because an approximation of a n to order 
n -m gj ves r j se ^ a difference equation for r n to order n~^ m+1 ^ via (J2T| . In turn, r n to 
order n _(m+1 ' ) produces an equation for a n to order n _ ^ m+1 ^ and so forth. This scheme 
breaks down for x ^ 0. 

For fixed i ^ 0, the numerics does not suggest an ansatz for the asymptotes. Most 
remarkably, the effect of 1 / persists for very large n, even for very small x 7^ 0, 
as illustrated in figure Also shown in this figure as dashed (dotted) lines are the 
approximations of a n (x) from the first three (two) terms of a Taylor-series in x around 
x = based on the iterative results for a n (x = 0), ()25|) and (J2Tj) . In principle, the 
Painleve IV, Equation (^7]) . provides a way to express a^ a \x) in terms of lower order 
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derivatives ai m \x) with m' < m, yet the results in figure [T] suggest that for sufficiently 
large n a finite Taylor series eventually deviates wildly from the correct a n (x). Note 
that both y/na n (x) and r n (x) are bounded in n for large n. 

Figure El shows the rescaled a n for fixed n and varying x. It resembles a Hermite 
polynomial because of its direct relation to P%(x,x)wq(x), (JTHJ) and (JTHJ. For the same 
reason a n (x) vanishes for nln(x) <C x 2 /2. 
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